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Abstract 

We study tt°tt correlations in the exclusive reaction pp — ► 4ir° at rest 
with complete reconstruction of the kinematics for each event. The inclusive 
distribution shows a dip at small invariant mass of the pion pair while a 
small enhancement in the double differential distribution is observed for small 
invariant masses of both pion pairs. Dynamical models with resonances in the 
final state are shown to be consistent with the data while the stochastic HBT 
mechanism is not supported by the present findings. 

1 Introduction 

Nucleon-antinucleon annihilation into multi-pion states offers the possibility of studying 
Bose-Einstein (BE) correlations under controlled conditions. While the study of BE corre- 
lations in inclusive distributions in NN annihilation in the conventional Hanbury-Brown- 
Twiss (HBT) framework M has a long history [Q-|4j, the use of exclusive distributions is 
relatively new. The first steps in this direction were made in (5|, ^| where the exclusive 
reactions pp — > 27r + 2-7r _ and pp — > 27r + 27r _ 7r° at rest were studied on the basis of min- 
imum bias CPLEAR data. For these data, the complete kinematical reconstruction of 
each event allows the direct determination of the square of the reaction amplitude and the 
dynamics of the BE correlations is not obscured by integration over spectators and the 
usage of conventional reference samples. 

The present paper extends the analysis of the exclusive distributions for the pp annihi- 
lation at rest to the case of the 4ir° final state on the basis of the Crystal Barrel (CB) data. 
An appealing property of the 47r° channel is that there are no p mesons, which strongly af- 
fect the final states with charged pions. Therefore the comparison of the 7r°7r° correlations 
with the charged pion correlations can clarify the origin of the signals observed. 

The HBT mechanism jlj, which is based on pion emission with stochastic phases 
over an extended emission volume, predicts an enhancement of like sign pion pairs at low 
relative momentum Q 2 . This kind of enhancement seen for tt + tt + and tt~tt~ pairs in many 
reactions is usually called a BE signal. For the first time in an annihilation reaction, the 
inclusive 2ir° pair distribution for pp — > Air shows no such enhancement contrary to the 



case of pp — ► 2-7r + 27r _ || where a weak enhancement was seen in the corresponding ir + ir + 
and 7r~7r~ distributions. In the double differential distribution, which is a more sensitive 
observable, a small enhancement is seen in the kinematical region where the 4-7r° system 
forms two pion pairs, both having small invariant masses. As we shall discuss, the strength 
and the shape of the signal for the differential correlations do not favor an interpretation 
in terms of the conventional HBT picture of BE correlations. This conclusion is based on 
simulations of resonance production in the 47r° final state which qualitatively explain the 
observed effects. 

The plan of the paper is as follows. In section ^ we describe the analysis of the data 
and the results for the single variable distributions. In section ||, the formalism of double 
differential distributions is recapitulated. Section || presents detailed model calculations 
for the dominant resonance mechanisms which are compared to the data. Partial projec- 
tions similar to the analysis performed earlier in the 27r + 27r~ case are shown as well and 
confronted with our dynamical model. Section [5| gives a summary and conclusions. 

2 Analysis of the 47r° data 

The Crystal Barrel experiment at the Low Energy Antiproton Ring (LEAR) at CERN 
was conceived for high statistics meson spectroscopy. The setup relevant for the present 
study used a cylindrical detector consisting of a hydrogen gas target, tracking chambers for 
charged particles (proportional and jet-drift-chambers), and an electromagnetic calorime- 
ter with 1380 Csl crystals of 15 radiation lengths each for the detection of photons. The 
calorimeter is surrounded by a 1.5 T solenoid magnet. 

2.1 Event selection 

The 47T° data analyzed in this study was taken by Crystal Barrel in November 1994 using 
antiprotons stopped in the hydrogen gas target at 12 bar with the "all-neutral trigger", 
i.e. only events with no signal in the proportional chambers were recorded. Offline, the 
data was subjected to the following selection criteria: 

• exactly 8 photon energy deposits (PED) in the calorimeter; 

• no PED adjacent to the beam pipe (risk of missing energy); 

• no p pileup condition; 

• no charged tracks in proportional chambers or jet-drift-chambers (eliminates backscat- 
ters, trigger inefficiency, and pair production); 

• no PED below 20 MeV(suppresses split-offs); 

• total momentum between and 200 MeV/c, and total energy mp p ±173 MeV. 

For this set of events, a kinematic fit is performed and the most likely photon combinations 
are chosen. Combinatorial background is determined by simulation to be of the order of 
18 percent. The photon energy resolution is cte/E = 2.4%/ tyE/GeV. 

The measured single pion inclusive momentum distribution is plotted in Fig.Q. Con- 
trary to the case of four charged pions 2tt + 2it~' ||], the shape of the tt° momentum distri- 
bution differs significantly from phase space due to the strong production of resonances in 
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Figure 1: The measured single-pion momentum distribution dN^/ d\p v \ for the Att° 
channel (thick histogram). The pure phase space distribution is shown with o and 
the single-pion momentum distribution for the 27r + 27r~ channel || with □. 



the mass region between 1.2 GeV and 1.7 GeV as discussed in [jUj]. The broad peak at 0.2 
GeV/c is due to pions recoiling against a heavy resonance like a2(1660) or 7^(1670), while 
the maximum around 0.65 GeV/c stems from the decay of lighter resonances like /2(1270) 
and maybe 7r(1300). The narrow peak at 0.85 GeV/c corresponds to the pp — ► tt° +r)_^ Sn o 
decay. We shall discuss the resonance production in detail in section |j. 



2.2 Correlation functions for inclusive distributions 

For the benefit of the reader, we summarize below the standard BE formalism according 
to ||. The single-particle inclusive density p\{p\) and the two-particle inclusive density 
P2(pi,P2) are related to the differential cross-sections by 

piiPi) = a "i^fm (1) 

. , _i da , . 

One of the definitions of pion pair correlations is based on the formula 

C{P1,P2) = P2(P1,P2) ~ Pl{Pl)Pl{P2) ■ (3) 

Alternatively the two-particle correlations can be described in terms of the ratio 

nl \ P2(P1,P2) , A s 

C(P1,P2) = —, r i 4 ) 

P0{Pl,P2) 

where po{pi,P2) is the two-particle distribution in the absence of correlations, with various 
prescriptions being used in the literature. A choice consistent with Eq.(|3|) is the product 
of the single-particle densities po(pi,P2) = Pi{pi)pi(P2)- 
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Averaging over angles and momenta gives a correlation function depending on one 
variable, the two-pion invariant mass M : 



c(M) = , (5) 

(pi-pi)(M) 

P2(M) = Js(m- sj^x + p 2 ?) p 2 (pi, P2 ) ^§|frj (6) 

(pi-pi)W = |* (M - ^(Pi +P2) 2 ) Pi(pi)pi(p 2 ) ^){2E 2 ) (7) 

The invariant mass M is uniquely related to the square of the momentum difference: 

{Pi-P2? = 4p 2 -M 2 = -Q 2 , (8) 

where p, is the pion mass and Q is the difference of the three-momenta of the two pions 
in their center-of-mass system (CMS), therefore the variables M 2 and Q 2 are equivalent. 

Because of the total energy-momentum conservation, the ratio C{M) is not a constant 
even if the distributions da/(d 3 p 1 /2Ei) and da/(d 3 p 1 /2Ei)(d 3 p 2 /2E 2 ) are determined by 
phase space alone (see HU). This effect becomes negligible for reactions at much higher 
energy, but it is important for annihilation at rest. 

2.3 Single-variable two-pion correlations 

In this subsection we present the single- variable two-pion correlation C(M) which has 
been frequently used in previous analyses. In order to isolate the correlation effects we 
compare the experimental density with a four-pion phase space distribution corrected for 
experimental cuts and efficiencies in the same way as the data. 

The data sample of 459803 events was used to calculate the two-particle distributions 
p 2 °(M) defined by Eq.@ for 7r° pairsQ, see Fig. ^a. The corresponding two-particle density 
from the phase space simulation is called p 2 (M). The simulation was produced with the 
GEANT software jnj, and subjected to the same cuts and selections as the real data. 
It comprises 956511 events. Figure 0b shows the ratio p 2 ° (M) / p 2 s (M) , for which the 



kinematical correlations discussed in section |2.2| cancel. An interesting feature of this 
result is that there is no familiar BE peak at small invariant mass of the pion pair, contrary 
to the case of pp — > 2ir + 2ir~ ||, which is shown for comparison in Fig. |2|b. The peak at 
M = 1.6 GeV is due to f 2 -> vrV. 

To study the correlation function C(M) in Eq.(||), the two-particle distribution for 
uncorrelated pion pairs has been calculated using the same event-mixing method as in ||. 
The experimental distribution {p\ ■ p\)(M) and the simulated distribution (p± ■ pi) PS (M) 
are plotted in figure ||, along with their ratio. Figure || shows the correlation functions 
C(M) for the 7r°7r° pairs in comparison with the phase space distribution. In order to 
remove the trivial M-dependence which arises from the energy-momentum conservation, 
equations (^-^), and the influence of the experimental cuts, the following double ratio has 
been calculated: 

C(M) _ p 2 (M) P ? S (M) 



C p S{M) (pi ■ pi){M) ■ (pi ■ pi) PS (M) 



^^Here and below all distributions for Air events contain multiple entries per event corresponding 
to all possible two-particle combinations. 
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Figure 2: Inclusive two-pion correlations vs. the effective mass squared M 2 of 
the pion pair in the 4.7T channel: (a) The experimental two-particle distribution 
P2{M), Eq.([]), in comparison with the phase space distribution p PS (M). The sim- 
ulated spectrum is normalized to the number of measured events. (b) The ratio 
p 2 (M)/p PS (M). The 27T+27T- data are from §. 
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Figure 3: a) The experimental distribution (pi • pi)(M) (•) and the phase 
space distribution (pi • pi) PS (M) (o) obtained by mixing events. b) The ratio 
(pi • Pi)(M)/(pi ■ pi) PS (M). The simulated spectrum is normalized to the number 
of measured events. The 2tt + 2'k~ data are from 0. 
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Figure 4: (a) The experimental correlation function C(M), Eq.([|), vs. the square 
of the effective mass M 2 of two pions in the 4ir° channel (•) in comparison with 
the corresponding phase space correlation function (o). (b) The experimental 
correlation function normalized to the phase space distribution, C(M)/C PS (M), 
Eq.(|). The 2n+2n- data are from §. 



The result is shown in FigJ^b. Clearly the inclusive ir° pair correlation function has a dip 
at small M. This is unexpected in the conventional HBT picture. It also contrasts to our 
findings for charged pion correlations, where the same correlation function shows a weak 
enhancement at small M ||, see Fig.[|b. As discussed in section |], the structures seen in 
FigfJ correspond to the resonances /2(1270), 02(1660), and 7T2(1670), 



3 Differential two-pion correlations 

So far we have presented the inclusive correlation function C{M)/C PS {M) where all kine- 
matical variables except the invariant mass M of one pion pair have been integrated out. 
To investigate possible correlation signals in more detail we turn to differential densities. 
Our approach follows the method used previously for the systems of charged pions f5[||]. 

The reaction pp — > 47r° at rest in the CB experiment proceeds from the 5-wave state 
J = ^ and the P-wave atomic states J PC = ++ , 1 ++ , 2 ++ . The corresponding pion 
distribution for the final state configuration {pj}, i = 1,2,3,4, has the form 

d<r({Pi}) ~ (^5|T(k,{ P J)| 2 + ^p|V fc T(k,{pJ)| 2 ) k ^ o d* 4 (p,Pi,P2,P3,P4) 

= |T(k, {Pj)|f^ d$4(p,Pl,P2,P3,P4) ■ (10) 

Here T(k, {Pj}) is the amplitude of the pp annihilation from the initial pp state with 
relative momentum k, ^$4(^,^1,^2)^3,^4) is the 4-body relativistic phase space, and the 
limit k — > implies the incoherent sumQ of the S and P-wave states with the corresponding 

2 The notation |T(k, {Pi})|? implies a sum over initial spin states J PC = h , ++ , 1 ++ , 2 ++ , 
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weights ws and iop. The four-vectors of the pions are pi = (Ei, Pj), and p = (2m p ,0) is 
the total four-momentum for pp annihilation at rest, m p being the proton mass. 




Figure 5: The two-pion subsystems a = (tti + 7r 2 ) and b = (n 3 + 7r 4 ) in the four pion 
final state. 

We introduce the two-pion subsystems a and b with four- momenta p a = (p\ +P2) and 
Pb = {P3 + P4) and invariant masses M a and M^, see hgure |5[ Given the invariant masses 
M a and Mj,, the double differential cross section is defined by integrating over the angles 
specifying the relative orientation of the momenta of the pions within the subsystems a 
and b and the relative orientation of the subsystems o and b (the corresponding solid angles 
are d£l\2, ^34, dQ ab ). For the reaction pp — > 4-/T we have 



da 



dM 2 a dMl 



w{^~s,M a ,M b ) J \T(k,{ Pi })\ 2 k ^ o dn ab dn 12 dn u . (11) 



Here the factor W(M,M a ,Mt,) is given by 



W(M,M a ,M b ) = ^ 



\ V Ml) 



4/48 H i-^l - (12) 



-Pal, = P(M,M a ,M b ) , (13) 

where P(M, M a , M b ) is the relative momentum of two particles with masses M a and M b 
and the total invariant mass M, in our case M = 2m p . Removing the phase space factor 
W(M, M a , M b ), we define the double-differential density: 

r , , „ , % 1 da 

g(M a ,M b ) ■- 



W(M, M a ,M b ) a-dMldMl 
~ I |T(k,{ Pi })|^ dQabdnnd^ (15) 

where a is the total cross section. 

An advantage of using the double differential density g(M a ,M b ) is that it does not 
contain kinematical dependences on the invariant masses of the two-pion pairs a and b. 
This means that for a constant T matrix the density g(M a ,M b ) does not depend on its 
argument M a and M b , contrary to p2(M a ) = da/dM a of Eq.(g) which unavoidably involves 
the phase-space dependence. 

The double-differential cross-section da/dM^dM 2 corresponding to the raw data is 
shown in figure [6|a, and the corresponding double differential density g{M a ,M b ) in fig- 
ure ||b. These raw data must be corrected for the acceptance of the detector, which is 



and all quantum numbers specifying the initial spin and the total angular momentum are sup- 
pressed. 
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Figure 6: The double-differential distributions vs. the invariant masses of two 7r° 
pairs (6 entries per one physical event): (a) the double-differential cross section 
da/dMldM%] (b) the double-differential density g(M a , M b ), Eq. The box 

plots (top) and the surface plots (bottom) represent the same data. Both the cross 
section and the density plots are not corrected for acceptance. 
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Figure 7: The acceptance simulation of the detector, calculated for events 
distributed uniformly in phase space: (a) the double-differential distribution 
da PS I 'dMldMl; (b) the double-differential density g PS (M a ,M b ). The box plots 
(top) and the surface plots (bottom) represent the same data. 
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shown in figure [j]. The differential density g PS (M a , M b ) is calculated using the four-pion 
phase space with the detector acceptance taken into account. The double differential den- 
sity corrected for the detector acceptance g(M a , Mj,)/g FS (M a , M&) is plotted in figure ^. 
Note that the double differential view shows a weak enhancement at low values of M 2 
which is absent in the inclusive projection of figure ||. 




Figure 8: The experimental double-differential density g(M a , M b )/ Q FS (M a , M b ) 
corrected for acceptance. Both plots represent the same data. 

For the further discussion in section ^| we construct the following partial projections 
of the two-particle density. The two-dimensional space (M a ,M b ) is divided into slices 
Mf < Mf, < M 2 +1 and the projections Qi(M a ) are defined by 

f M ? + i g(M a ,M b ) a 

e ^ Ma) = l Mf 7HM-m dMb (16) 

The partial projections Qi(M a ) are plotted in Fig. ||. Contrary to the inclusive distribution 
figure |I]b where a dip was observed, we see a peak in the first two slices. The slices in the 
mass interval 0.3 GeV 2 < M 2 < 1.2 GeV 2 are more or less flat. The dip seen in Fig. |]b 
results from integration over the higher mass region of M 2 > 1.2 GeV 2 . In comparison to 
the corresponding projections for the annihilation into 27r + 27r _ , the observed enhancement 
at small Q 2 in the irir case is weak even for the lowest interval 0.0 < Q 2 < O.lGeV 2 . 



4 Discussion 

4.1 Resonance effects 

Resonances are strongly produced in the reaction pp — » 4tt°, see figures |6] where the 
signals of fa are clearly seen in the invariant mass projections. The total spin S and the 
relative angular momentum L of the pp system annihilating into 47r° are constrained by 
the conservation of C-parity: C = (— \) L+S = 1. Therefore, the annihilation into Air 
can occur in the case of the 5-wave annihilation (L = S = 0) from the state J PC = 
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Figure 9: The acceptance-corrected projections Qi(M a ) of the differential density 
g(M a , Mb) I g PS (M a , Mb) for different intervals of the invariant mass of the other 7r°7r° 
pair. The relative momentum Q, Eq.(H), is used instead of the invariant mass M 
for the sake of convenience. 
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only, and in the case of the P-wave annihilation (L = S = 1) from the states J PC = 
0++,l++,2++. 

Following the Crystal Barrel analysis []TI| we focus our attention on the most prominent 
mechanisms seen in the Air channel. For the dense gaseous target at 12 bar, the fractions 



of the total S- wave and P-wave annihilations are approximately equal to each other [15]. 
For the sake of simplicity, we focuse our attention on one P-wave channel J = 2 ++ , 
which has the largest statistical weight, and the S- wave channel J PC = h . We do not 
attempt to make a global fit of the data in the Air channel, which would be far beyond the 
scope of this paper, but rather to explore how the most probable annihilation mechanisms 
manifest themselves in the observed BE correlations. The following mechanisms have been 
taken into consideration: 

pp(J PC = (n + ) -► a 2 (1660)+^° (17) 

^/ 2 (1270)+vr° 
^2tt° 

pp(J PC = 2++) -» / 2 (1270) + 2tt (18) 

^2tt° 

pp(J PC = 2 ++ ) -► a 2 (1660)+^° (19) 



/ 2 (1270)+vr 



o 



27T 



pp(J PC = 2 + +) vr 2 (1670)+^° (20) 

^/ 2 (1270)+vr° 



2vr° 



pp(J PC = 2++) 7r(1300)+ 7T° (21) 

a + vr 



2tt° 



pp(J PC, = 2++) - a + fo (22) 

^4vr° 

pp(J PC = 2 ++ ) -» a + a (23) 

p p (J PC = 0++) a + a (24) 

^4tt° 

pp(J PC = 0++) v + fo (25) 



^4vr° 

The details of the calculations are given in Appendix [A|. 

As a first step, we inspect the double-differential densities corresponding to each indi- 
vidual mechanism in order to demonstrate the effects expected if one of these mechanisms 
would be dominant. In reality, a coherent superposition of several mechanisms is required 
to describe the data pX] , therefore, as a second step, we investigate the double-differential 
distributions corresponding to such superpositions. Our goal is to find out which mech- 
anisms are required to describe the global features of the differential mass distributions. 
Once these features are established we also get a prediction for the double differential dis- 
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Figure 10: The double-differential density g(M a , Mb) calculated for different reac- 



tion mechanisms: (a) pp(J PC = + ) 
2++) - / a (1270) + 2vr , (c) p P (J PC 
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Figure 11: The double-differential density g{M a ,M b ) calculated for a sum of the 
mechanisms (]T^-f2D|) as discussed in the text. 



tributions near the origin from these ordinary resonance mechanisms without introducing 
any HBT-like stochastic correlations. 

The double-differential densities calculated for the mechanisms (17,18,19,^) are shown 
in Figs.[l(](a,b,c,d). 

By varying the relative weights and phases of the individual mechanisms we have found 
that the superpositions of the three mechanisms (15,19, and 20) are sufficient to obtain a 
qualitative agreement with the data: 



ClT /2(1270)+2vrO +C2d l(t>2 T a2 (i 66 o) +7r O + C 3 e i<p3 T^ lim ) 



i<t>3 



(26) 



In fact, the main features of the data can be reproduced even with the two mechanisms 
fll9|) and (pc|) (note that they include the contribution from the lighter meson /2(1270), 
see Figs.[l2|(b,c)). Figure 11 shows an example corresponding to a superposition of the 

3vr/2 (all 



mechanisms ( |19|) and (|20[) with c\ = 0, C2 = 0.45, C3 = 0.91, (f>2 = 0, 03 
individual amplitudes are normalized to the same total yield of 4vr°). We are not aiming 
at an optimal fit to the data in the present context, but Fig. [ll] illustrates that the observed 
resonance peaks are consistent with a moderate enhancement near the origin. Compare 
this figure with the experimental distribution of figure 0. 



5 Conclusion 

For the annihilation reaction pp — ► 4n° at rest, the BE correlations between two neutral 
pions have been studied for the first time with full event reconstruction. The inclusive 
27T° correlation shows an unexpected dip at small relative momentum, in contrast to 
the inclusive pair correlation in the charged channel pp — > 2tt + 2tt~ which shows a weak 
enhancement ||. 

For the more sensitive and less model-dependent double-differential distribution, a 
weak peaking is observed near the origin (small invariant masses of the two neutral pion 
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pairs). This is similar to though less pronounced than in the 2tt + 2'k~ channel. 

The data of the single- and double-differential 2ir° correlations were compared to model 
calculations. The role of the p meson dominating the 2it + 2it~ dynamics is replaced by 
the d2, 7T2, and fa resonances in the 4-7T case. Resonance production explains the fairly 
different spectra of the two reactions. It is found that with an adequate choice of reaction 
amplitudes the dynamical model describes the main features of the single- and double- 
differential distributions, including a slight enhancement near the origin. 

From this result, together with the absence of a correlation signal in the inclusive 
2n° correlation, we conclude that our analysis does not favour the interpretation of pion 
correlation signals in terms of an HBT-type model with stochastic pion emission phases. 

6 Aknowledgements 

The authors are very grateful to the Crystal Barrel Collaboration for making available the 
data. In this connection we would like to thank in particular Ulrich Wiedner and Crtomir 
Zupancic. We had stimulating discussions on aspects of pion correlations with Eberhard 
Klempt and Ulrike Thoma. One of us (MPL) would like to thank the TRIUMF theory 
group for hospitality during completion of this paper. 



A Appendix 

The amplitudes corresponding to the resonance mechanisms (|]||ijj2||||[i]].||) are given 



by the diagrams in figures 12. The calculations are straightforward, so only some ma- 
jor elements are outlined here. The pp vertices are taken with a minimum number of 
derivatives. The vertex pp(J PC = 2 ++ ) — » a2vr has the form 

Vp p{J pc =2++) ^ a27T = B, u g^T lS e a ^ s {pM Pa2 ) p (27) 

where is the polarization tensor of the pp( J PC = 2 ++ ) state and T^s is the polarization 
tensor of the a-i, and p^ and Pa2 are the momenta of the corresponding particles. The 
0-2 fa^ vertex has a similar structure, and the pp(J PC = 2 ++ ) — * 1T2TT is a trivial 
5-wave vertex. The fawn vertex has the form 

V)^ = 9f 2 -!r7rT^q tl q u Ff 27nr (ml^) (28) 

where j/ 2 „ is the coupling constant, T^ v is the polarization tensor of the fa, q = p\ — P 2, 
p n are the pion momenta, and Ff 27r7T (m^ n ) is the form factor depending on the invariant 
mass of the tttt system. 

The propagators corresponding to the tensor particles are given by the formula 

i(n^'n w ' +W^'W V ') - ±(IT"TP V ) 
G{P) ~ pt-ml-Mb?) (29) 

IT"' = ^-/p7p 2 (30) 

where mo is the bare mass and M{p 2 ) is the mass operator. The mass operator for the fa 
meson corresponding to the tttt loop is defined by the dispersion integral: 

= i r n*l ds , (31) 

2vr J 4m 2 s -s' 

T(s) = g^Wflg (32) 

607TS 
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where k = k(s) = \J s/A — m 2 is the relative momentum in the tttt system and the dipole 
form factor Ff 27Tn (s) = (1 + k(s) 2 /u 2 )^ 2 is used with v = 1 GeV. The parameters mo and 
gf 27T7T are defined by the mass and the width of the f2 meson, with the other decay channels 
being approximated by a constant width. The propagators of 02 and 7T2 are constructed 



in a similar way. The a block in Fig. 12 a denotes the full Green function of the tttt system 



in the scalar- isoscalar channel; it was taken from the coupled channel model [16|. The 
components above are combined to construct the amplitudes for the resonance channels 
concerned; since the procedure is straightforward, we skip the resulting formulas involving 
lengthy tensor structures. 
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